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Abstract 

The E* — E c and E£ — E;, hyperfine mass splittings are computed in the Skyrme 
model. The hyperfine splittings are suppressed by both 1/N C and by 1/mg, where N c is 
the number of colors and ttiq is the mass of the heavy quark. The E c , E*, E£, and A;, 



^ . masses are predicted in terms of the known values of the A c , D, D* , B and B* masses. 
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Baryons containing a heavy quark can be interpreted as bound states of mesons con- 
taining a heavy quark and chiral solitons of the SU(2)l X SU(2)r nonlinear sigma model. 
This approach was originally proposed by Callan and Klebanov to treat baryons contain- 
ing an s quark as soliton-iiT-meson bound states |l| . Early attempts were made to extend 
this method to baryons containing a heavy quark Q, such as the c or b quark MP]- A 
correct treatment of heavy baryons as soliton-heavy-meson bound states must incorporate 
the consequences of heavy quark symmetry ||. In the heavy quark limit, the D and D* 
(B and B*) mesons are degenerate, so the effective chiral theory must include both the 
pseudoscalar (Pq) and vector meson (Pq) fields in order to respect the heavy quark sym- 
metry. Properties of the heavy baryon bound states have been calculated to leading order 
in 1/mq @0- The Skyrme model successfully predicts the existence of Ag, Eg and Eg 
bound states and leads to mass relations which are well-satisfied by the measured A c ^ 
and £ Ci b masses [|J. In the heavy quark limit, the Eq and Eg are a degenerate multiplet 
with isospin one, spin of the light degrees of freedom one, and total spin 1/2 and 3/2, 
respectively. The energy splitting between these states is a 1/mq effect. In this letter, the 
Eg — Eg mass difference is computed. 

It is instructive to first study the hyperfine splittings in the constituent quark model. 
The Eg — Eq and Pq — Pq mass differences are due to the hyperfine interaction generated 
by one-gluon exchange between constituent quarks [BJ, 



number of colors, and the sum is over all pairs. The constant k depends on the wave 
function of the hadron, and may be different for mesons and baryons. 7i has an overall 
factor of 1/N C , since the gluon coupling constant is of order 1/V iV c . The pseudoscalar 
and vector mesons Pq and Pq are bound states of a heavy quark Q and a light antiquark 
q in a color singlet state, with spin zero and one, respectively. The hyperfine interaction 
Eq. (TJ) in the meson sector can be written in the form 




(1) 



where T^, Si and rrii are the color generator, spin and mass of the z th quark, iV c is the 
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for the SU(N C ) representation R. This gives the hyperfine splitting 
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where the irreducible representations of SU(N C ) are denoted by Young tableaux. 

The hyperfine interaction Eq. ([!]) in the baryon sector can be written in the form 
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where the first term is the interaction of the heavy quark with the light quarks, and the 
second term is summed over all light quark pairs. The last Casimir is that of the two-index 
antisymmetric tensor, because any two quarks in a baryon are antisymmetric in color. The 
spin of the light degrees of freedom is the spin of the light quarks, Sg = Yli^i 1 Si, so that 

n is 

Sq ■ Si 



n = 



K 

2W r 



2c(n)- c (f]; 



_m Q m q 



(6) 



where Tig is the interaction energy of the light degrees of freedom. Tig is of order one, since 
all the explicit N c dependence has been factored out. The Eq and Eg have light degrees 
of freedom in identical quantum states so that Tin cancels in the mass difference. Using 
2S Q ■ S e = (S Q + Sgf -S 2 Q - Sj, with sq = 1/2, Si = 1, and s Q + s e = 1/2 for the Eg 
and sq + si = 3/2 for the Eg, gives 

„ 3k' 
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The Ag state is a spin-1/2 baryon with the spin of the light degrees of freedom sg = 0, so 
the first term in Eq. (^|) does not contribute to the Ag mass. The first term in Eq. (|B|) also 
cancels out of the average mass of the Eg, Eg multiplet (Eg) avg = -|(Eg + 2Eg). The 
(Eg) av g — Ag mass difference is 



(EQ) av g - Ag 
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where ATig is the difference in interaction energies of the light degrees of freedom in the 
Ag and Eg, Eg. 

The Casimirs are trivial to compute using standard techniques of group theory M, 



c (i) = o, c (n: 



_ N 2 -l 
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(iV c -2)(iV c + l) 



(9) 



The color factor for the mesons is 



_ N 2 -1 
c(D) + c(D)-c(l) = -^ >N C (N c ^oo), (10) 
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and for the baryons is 

2c(n)- C (g) = ^2+l-l (JVc 



Oo). (11) 



Thus, in the large iV c limit, the hyperfine splittings are 

K 



Pq-Pq 



2mQm q 
k' 

(Sg)avg - Aq = ^rAHt. 

The Pq — Pq mass difference is of order 1/toq, the Eg — Eq mass difference is of order 
(l/N c )(l/m,Q), and the (EQ) avg — Aq mass difference is of order 1/N C . 

The Eg — Eq and (Eg) avg — Aq mass differences can be computed in the Skyrme 
model, in which baryons containing a heavy quark are treated as soliton-heavy-meson 
bound states. The formulas that are needed from earlier work are summarized here; the 
details can be found in Refs. 0[[§. The Pq and Pq states are combined into a heavy 
meson field 

Ha = [P a *,7 M - Pal,] , (13) 

where a = 1, 2 denotes the Qu and Qd states, and v M is the heavy quark four- velocity. The 
transformation rule for the H field under SU{2)q heavy quark spin symmetry is 

H a -> SH a , (14) 

where S G SU(2)q, and the transformation rule under chiral SU{2)l X SU(2)r is 

H a ^{HB)) a , (15) 

where we use the primed basis for the H fields defined in Ref. The Goldstone boson 
fields have the SU{2)l x SU(2)r transformation law 

E(x) -> L E(x) P f . (16) 
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The E field can be written in terms of the pion fields as 



= e 2iM / f , (17) 



where M is the Goldstone-boson matrix 



M = 



TV°/V2 7T + 

7T" -n°/^/2 



(18) 



and the pion decay constant / ps 132 MeV. The soliton solution of the SU(2) L x SU(2) R 
chiral lagrangian is 

E = AE A-\ (19) 

where 

E = exp[iF(|f|) x -f], (20) 

and A e SU(2) is the collective coordinate associated with isospin transformations of the 
soliton solution £q- The radial shape function F(r) satisfies F(0) = — 7r and F(oo) = for 
a soliton with baryon number one. In the quantum theory, baryons have wavefunctions 
that are functions of the matrix A e SU{2). The wavefunctions are 



IpRam(A) = (-l) R+m Vdh^R Dl^UA), (21) 

for a state with isospin / = R, Is = a, J = R, and J3 = m. The matrices Z)^^ are 
the representation matrices for SU(2), and we have normalized the measure on the SU{2) 
group so that 

dA = 1. (22) 



Isi 



<SU(2) 

In QCD, the only soliton states in the theory are those with 21 = 2J = odd. 

In the large N c limit, the soliton is very heavy, and the semiclassical approximation 
is valid. Each time derivative is suppressed by a factor of 1/N C , and can be neglected to 
leading order in N c . The interaction hamiltonian to leading order in iV c and in the spatial 
derivative expansion is 

H T = -^ [d s x TrM7^ 75 [Et<9 M E], 

\ (23) 
w - y / d 3 x TrM^^sptdjE], 



with E given by Eq. ([19]). In the limit that the H field is very heavy, the interaction 
energy is determined by the value of Eq. (^3|) with the H field at the origin Q . Using the 
expansion 

F(r) = F(0) + rF'(O) + . . . = -tt + rF'(O) + . . . , (24) 
in Eq. ( |23| ) gives the interaction hamiltonian 

Hi = gF ^ / rf 3 x Tr HH 1 j l5 r k TrA-r^V, 

4 7 (25) 

= ^F'(0) 4 S{ H Tr,WV, 

where Ijj and Sm are the isospin and spin of the light degrees of freedom of H. The Eq, 
Eq multiplet has 7 = 1 and se = 1. The soliton states have I = si = 1/2, 3/2, . . ., and 
the light degrees of freedom of H have I = si = 1/2. There are two states in the soliton- 
meson spectrum with I = s# = 1, obtained by tensoring the light degrees of freedom of 
the H particle with either the iV (I = S£ = 1/2) or the A (J = s$ = 3/2) soliton states. 
These states will be denoted by |1 1; N) and |1 1; A), where the first label refers to isospin 
and the second to the spin of the light degrees of freedom. The interaction hamiltonian 
Eq. (|25|) in the |1 1;N, A) basis can be evaluated explicitly in terms of Clebsch-Gordan 
coefficients [|J, 

9 F'(0) ( 1 4v/2\ 
Hl = 6~U/5 5 J' (26) 

The two |1 1) states which are the eigenstates of Hi in Eq. (EBI) are 



|i i» = v| |ii;JV) + yf in; ^ 
11% = a/1 lii;Jv>-\/| |11;A) 



(27) 



with energies — 3(7i ?/ (0)/2 and gF'(0)/2, respectively. The state |1 1)) when combined 
with the heavy quark in H is the Eq, Eq multiplet, and |1 combined with the heavy 
quark is an unbound exotic multiplet. The subscript on the | )) states is the eigenvalue of 
K = I + Si; the significance of K is discussed in detail in Ref. 0]. The Aq has I = si = 0. 
The only state in the soliton-meson spectrum with the correct quantum numbers is in 
the nucleon sector, |0 0; N). The interaction energy of this state due to the hamiltonian 
Eq. (H) is -3gF'(0)/2 @. The state |0 0; N) has K = 0, and will also be denoted by 
|0 0» . 
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The Eg, Eq and Aq states are obtained by combining the spin of the light degrees 
of freedom with the spin of the heavy quark. Since spin states with different values of S3 
all have the same mass, only the Eq state with S3 = 3/2, the Eq state with s 3 = 1/2, and 
the Aq state with S3 = 1/2 are needed to determine the energies, 

|Sq, §>=|1H»oIT>o 

|S Q , \) = ^1 |1 1 1» U>g - |1 1 0» |t)g (28) 
|A Q , i>= |0 0» |t>g. 

In Eq. (|28|) , the state | )q denotes the spin state of the heavy quark, and \I sg m)) 
represents the state of the light degrees of freedom \I sg)) with = m. The coefficients 
in Eq. (|28|) are the Clebsch-Gordan coefficients for 1 <g> 1/2 -»• 3/2, 1 ® 1/2 -> 1/2, and 



(8) 1/2 -v 1/2, respectively. 

The leading term in the 1/N C and spatial derivative expansion that splits the N, A 
degeneracy is the kinetic term of the soliton, which has two time derivatives and no space 
derivatives, and is of order 1/N C . Including this term shifts the energy of the |1 1;A) 
states, and changes the interaction hamiltonian Eq. (I2B) to 



_ gF^f 1 4>/2\ , (0 1 
Hl ~ 6 1 4^/2 5 J + U AM'' (29) 



where AM is the A — N mass difference, and is of order 1/N e . The eigenstate of Eq. (|29| ) 
that (combined with the heavy quark states) produces the Eq, Eq multiplet and reduces 



to |1 1)) in the limit AM -> is 



|1 l)) e = a |1 1; N) + b |1 1; A) = |1 1)) + e |1 , (30) 



where 



1 12 , [2 fl AM . . 

3 + vr- " = V5-vr- e= 3V2 9 F'(o)' (31) 

to first order in AM. The Eq, Eq states to order 1/N C are given by Eq. (|2^ ) with |1 1)) 
replaced by |1 l)) e . Since the Aq state |0 0)) only contains the soliton in the nucleon 
sector, its interaction energy — 3gF' (0)/2 is unaffected by the A — N mass difference. 



The soliton-meson states |1 1; N) and |1 1; A) which appear in Eq. ( |30|) can be written 
explicitly in terms of the soliton states and the spin of the light degrees of freedom of H, 

|iii;^) = ll>^it),, 

l 11 0;iV) = v /||i> 7V |I), + v /||-i>^|T),, 

|iii;A) = v /|||> A u),-y / |||> A |T), ) (32) 

|iiO;A) = v /||i) A ||),- v /||-i> A |T),, 

\ooo',m = ^l\i) N \i)t-^l\-i) N n)t, 

where | ) e is the spin state of the light degrees of freedom in the heavy meson, \m) N A is 
the soliton state in the N or A sector with = m, and \I S£ m; N, A) is the bound state 
| J si; N, A) with S£3 = m. The coefficients in Eq. ( |3~2"D are the Clebsch-Gordan coefficients 
for 1/2 <g> 1/2 -> 1, 3/2 ® 1/2 -»• 1 and 1/2 ® 1/2 -> for the |1 1; JV), |1 1; A) and |0 0; iV) 
states, respectively. Finally, the tensor product of the light degrees of freedom and the 
heavy quark in H can be reexpressed in terms of P and P* mesons, 

U)*U>q = |*5.-i>. 

Wi\i) Q = \l\ \ p Q) + \fl I p q>°>> (33) 

i^iT) Q = -/|iPQ) + y / ||p^°>» 

where \Pq, to) is the Pg meson with S3 = m. Combining Eqs. (|2~8|)-(|3~3|), gives 

1% l> = -vf & II>aI*q)-^ II>aI^ 1 > + « IDivl^ 1 ) 

+ /! 6 ll> A l*Q.0>, 
|Sq, |> = yf a |i>„|P > + ^& ||> A |^,-l> + y|& I-DaI^ 1 ) ^ 

|Aq, \) = \ \\) n \Pq^)-\ \\) n \Pq)-\[\ \~\) n \Pq^)- 
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The leading operator in the derivative expansion which breaks the heavy quark sym- 
metry is 

& = — Tr -Ha^Haau (35) 
m Q 

which yields a Pq — Pq mass difference — 8X2/ tuq. The Eg, Eg, Ag masses including 
Eq. (|3^) can be read off directly from the expressions for the states in Eq. (|34"D, 

^ = -& + a ^ + b ^ + ^ PQ+ (^ + b ^ Pi , (36, 

where the coefficient a 2 is the probability that Eg contains a nucleon, 3o 2 /8 is the proba- 
bility that Eg contains a Pq, etc. The mass differences are 



3 [2a 2 - b 2 } . 

Zq-z q = ^- 8 l(p*-p Q ), 

(Eg) avg - A Q = b 2 (A - N) , 
using a 2 + b 2 = 1. Substituting the values of a and b from Eq. (|3l| ) gives 

v * v _ (A-AQ (Pe-P ) 

(E Q ) avg - Ag = I (A - TV) . 



(37) 



(38) 



Note that the mass splittings have the correct dependence on mg and N c . The Pq — Pq 
mass difference is of order 1/mg, and the A — N mass difference is of order 1/N C . This 
implies that the Eg — Eg mass difference is of order l/(mgiV c ) and the (Eg) avg — Ag 
mass difference is of order 1/N C , as expected from consideration of the constituent quark 
model. 

The value of gF'(0) is obtained from the A c mass in Eq. (|36|) to be 418 MeV, so that 



E* — E c = 25 MeV, E£ - E 6 = 8 MeV, 

(39) 

(E:) avg - A c = (E*) avg - A, = 195 MeV. 



This gives the predictions 

(input) A c = 2285 MeV A b = 5625 MeV, 

E c = 2464 MeV, E 6 = 5816 MeV, (40) 
E* = 2489 MeV, E£ = 5824 MeV, 

to be compared with the experimental values E c = 2453 MeV, and A& = 5641 ± 40 MeV. 
The predictions are in good experiment with experiment. Note that the typical magnitude 
of the E c — A c splitting is of order the A — N splitting, which is ~ 300 MeV. The Skyrme 
model predicts the value of this splitting to within « 10 MeV. The masses of the charm and 
bottom baryons have also been computed in the constituent quark model fTj. The quark 
model model prediction of the Eg — Eg splitting is about twice as big as the Skyrme model 
prediction. 

The hyperfine splittings in Eq. ([55p have been computed to leading order in the 1/N C 
expansion, and to leading order in the spatial derivative expansion, so they are not true 
large N c predictions. For example, the operator 

£ 2 = ^TrH^ 5 H^d^, (41) 
m Q 

violates the heavy quark symmetry, and destroys the equality between the PqPqTt and 
PqPqtt coupling constants. This operator will produce a Eg — Eg mass difference pro- 
portional to (A 7r /mg)(A — N). The operator £ 2 has one space derivative, whereas L\ in 
Eq. fl35| ) has no space derivative. If the derivative expansion is valid, the contribution of 
£2 to the Eg — Eg mass difference should be smaller than the L\ contribution that we 
have computed. 

Some general features of the results in this paper can be easily understood on the 
basis of K invariance. In the large N c limit, baryon states can be classified by the quantum 
number K, and the Ag and Eg, Eg states all have K = 0. At order 1 /N c the Eg, Eg states 
are a superposition of a K = state with amplitude N® and a K = 1 state with amplitude 
1/N C . The Ag state remains a pure K = state because there is no K = 1 state with 
the correct quantum numbers to mix with the Ag. The Eg — Eg mass difference violates 
the heavy quark symmetry, and therefore must come from an operator in the effective 
hamiltonian that involves Sq, the spin of the heavy quark. By rotational invariance, Sq 
must be dotted into the spin of the light degrees of freedom of either H or the soliton. 
(The orbital angular momentum contributions are suppressed by 1/N C or by an additional 
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1/rriQ.) The operator that breaks the heavy quark symmetry must therefore transform as 
a K = 1 irreducible tensor operator, since K = I + Se- A K = 1 tensor operator cannot 
have a matrix element between two K = states. Thus the operator that breaks the heavy 
quark symmetry does not affect the Aq energy, and only contributes to the Eq, Eq energy 
at order 1/N C . The 1/mg terms such as Eq. (^5j) are of leading order in 1/N C because they 
contain no time derivatives. The 1/N C K- violation arises from the A — N mass difference 
in the soliton sector. There also are terms that simultaneously violate the heavy quark 
symmetry and the X-symmetry, and are of order 1/(itiqN c ) because they contain one time 
derivative, e.g. 

C 3 = —e^BTrH^^H^d^d^d^ (d 2 E+E - E^E) , (42) 

or 

C 4 = — Tri5V 75 H 7^75 E f (v ■ d) E (d 2 E f E - E f a 2 E) . (43) 

m Q 

Both terms start at one time derivative and four space derivatives, and are expected to 
be small. (£4 naively starts at two space derivatives, but Y^dH/dt vanishes at the origin 
where the H particle is bound, which brings in an extra suppression factor.) 
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